\(\int \frac {d+e x^2}{(a+c x^4)^2} \, dx\) [146]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [C] (verified)
   Fricas [B] (verification not implemented)
   Sympy [A] (verification not implemented)
   Maxima [A] (verification not implemented)
   Giac [A] (verification not implemented)
   Mupad [B] (verification not implemented)

Optimal result

Integrand size = 17, antiderivative size = 275 \[ \int \frac {d+e x^2}{\left (a+c x^4\right )^2} \, dx=\frac {x \left (d+e x^2\right )}{4 a \left (a+c x^4\right )}-\frac {\left (3 \sqrt {c} d+\sqrt {a} e\right ) \arctan \left (1-\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt {2} a^{7/4} c^{3/4}}+\frac {\left (3 \sqrt {c} d+\sqrt {a} e\right ) \arctan \left (1+\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt {2} a^{7/4} c^{3/4}}-\frac {\left (3 \sqrt {c} d-\sqrt {a} e\right ) \log \left (\sqrt {a}-\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{16 \sqrt {2} a^{7/4} c^{3/4}}+\frac {\left (3 \sqrt {c} d-\sqrt {a} e\right ) \log \left (\sqrt {a}+\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{16 \sqrt {2} a^{7/4} c^{3/4}} \]

[Out]

1/4*x*(e*x^2+d)/a/(c*x^4+a)-1/32*ln(-a^(1/4)*c^(1/4)*x*2^(1/2)+a^(1/2)+x^2*c^(1/2))*(-e*a^(1/2)+3*d*c^(1/2))/a
^(7/4)/c^(3/4)*2^(1/2)+1/32*ln(a^(1/4)*c^(1/4)*x*2^(1/2)+a^(1/2)+x^2*c^(1/2))*(-e*a^(1/2)+3*d*c^(1/2))/a^(7/4)
/c^(3/4)*2^(1/2)+1/16*arctan(-1+c^(1/4)*x*2^(1/2)/a^(1/4))*(e*a^(1/2)+3*d*c^(1/2))/a^(7/4)/c^(3/4)*2^(1/2)+1/1
6*arctan(1+c^(1/4)*x*2^(1/2)/a^(1/4))*(e*a^(1/2)+3*d*c^(1/2))/a^(7/4)/c^(3/4)*2^(1/2)

Rubi [A] (verified)

Time = 0.13 (sec) , antiderivative size = 275, normalized size of antiderivative = 1.00, number of steps used = 10, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.412, Rules used = {1193, 1182, 1176, 631, 210, 1179, 642} \[ \int \frac {d+e x^2}{\left (a+c x^4\right )^2} \, dx=-\frac {\arctan \left (1-\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right ) \left (\sqrt {a} e+3 \sqrt {c} d\right )}{8 \sqrt {2} a^{7/4} c^{3/4}}+\frac {\arctan \left (\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}+1\right ) \left (\sqrt {a} e+3 \sqrt {c} d\right )}{8 \sqrt {2} a^{7/4} c^{3/4}}-\frac {\left (3 \sqrt {c} d-\sqrt {a} e\right ) \log \left (-\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {a}+\sqrt {c} x^2\right )}{16 \sqrt {2} a^{7/4} c^{3/4}}+\frac {\left (3 \sqrt {c} d-\sqrt {a} e\right ) \log \left (\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {a}+\sqrt {c} x^2\right )}{16 \sqrt {2} a^{7/4} c^{3/4}}+\frac {x \left (d+e x^2\right )}{4 a \left (a+c x^4\right )} \]

[In]

Int[(d + e*x^2)/(a + c*x^4)^2,x]

[Out]

(x*(d + e*x^2))/(4*a*(a + c*x^4)) - ((3*Sqrt[c]*d + Sqrt[a]*e)*ArcTan[1 - (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/(8*Sqr
t[2]*a^(7/4)*c^(3/4)) + ((3*Sqrt[c]*d + Sqrt[a]*e)*ArcTan[1 + (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/(8*Sqrt[2]*a^(7/4)
*c^(3/4)) - ((3*Sqrt[c]*d - Sqrt[a]*e)*Log[Sqrt[a] - Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2])/(16*Sqrt[2]*a^(
7/4)*c^(3/4)) + ((3*Sqrt[c]*d - Sqrt[a]*e)*Log[Sqrt[a] + Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2])/(16*Sqrt[2]
*a^(7/4)*c^(3/4))

Rule 210

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(-(Rt[-a, 2]*Rt[-b, 2])^(-1))*ArcTan[Rt[-b, 2]*(x/Rt[-a, 2])
], x] /; FreeQ[{a, b}, x] && PosQ[a/b] && (LtQ[a, 0] || LtQ[b, 0])

Rule 631

Int[((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> With[{q = 1 - 4*Simplify[a*(c/b^2)]}, Dist[-2/b, Sub
st[Int[1/(q - x^2), x], x, 1 + 2*c*(x/b)], x] /; RationalQ[q] && (EqQ[q^2, 1] ||  !RationalQ[b^2 - 4*a*c])] /;
 FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]

Rule 642

Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Simp[d*(Log[RemoveContent[a + b*x +
c*x^2, x]]/b), x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[2*c*d - b*e, 0]

Rule 1176

Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[2*(d/e), 2]}, Dist[e/(2*c), Int[1/S
imp[d/e + q*x + x^2, x], x], x] + Dist[e/(2*c), Int[1/Simp[d/e - q*x + x^2, x], x], x]] /; FreeQ[{a, c, d, e},
 x] && EqQ[c*d^2 - a*e^2, 0] && PosQ[d*e]

Rule 1179

Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[-2*(d/e), 2]}, Dist[e/(2*c*q), Int[
(q - 2*x)/Simp[d/e + q*x - x^2, x], x], x] + Dist[e/(2*c*q), Int[(q + 2*x)/Simp[d/e - q*x - x^2, x], x], x]] /
; FreeQ[{a, c, d, e}, x] && EqQ[c*d^2 - a*e^2, 0] && NegQ[d*e]

Rule 1182

Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[a*c, 2]}, Dist[(d*q + a*e)/(2*a*c),
 Int[(q + c*x^2)/(a + c*x^4), x], x] + Dist[(d*q - a*e)/(2*a*c), Int[(q - c*x^2)/(a + c*x^4), x], x]] /; FreeQ
[{a, c, d, e}, x] && NeQ[c*d^2 + a*e^2, 0] && NeQ[c*d^2 - a*e^2, 0] && NegQ[(-a)*c]

Rule 1193

Int[((d_) + (e_.)*(x_)^2)*((a_) + (c_.)*(x_)^4)^(p_), x_Symbol] :> Simp[(-x)*(d + e*x^2)*((a + c*x^4)^(p + 1)/
(4*a*(p + 1))), x] + Dist[1/(4*a*(p + 1)), Int[Simp[d*(4*p + 5) + e*(4*p + 7)*x^2, x]*(a + c*x^4)^(p + 1), x],
 x] /; FreeQ[{a, c, d, e}, x] && NeQ[c*d^2 + a*e^2, 0] && LtQ[p, -1] && IntegerQ[2*p]

Rubi steps \begin{align*} \text {integral}& = \frac {x \left (d+e x^2\right )}{4 a \left (a+c x^4\right )}-\frac {\int \frac {-3 d-e x^2}{a+c x^4} \, dx}{4 a} \\ & = \frac {x \left (d+e x^2\right )}{4 a \left (a+c x^4\right )}+\frac {\left (\frac {3 \sqrt {c} d}{\sqrt {a}}-e\right ) \int \frac {\sqrt {a} \sqrt {c}-c x^2}{a+c x^4} \, dx}{8 a c}+\frac {\left (\frac {3 \sqrt {c} d}{\sqrt {a}}+e\right ) \int \frac {\sqrt {a} \sqrt {c}+c x^2}{a+c x^4} \, dx}{8 a c} \\ & = \frac {x \left (d+e x^2\right )}{4 a \left (a+c x^4\right )}+\frac {\left (\frac {3 \sqrt {c} d}{\sqrt {a}}+e\right ) \int \frac {1}{\frac {\sqrt {a}}{\sqrt {c}}-\frac {\sqrt {2} \sqrt [4]{a} x}{\sqrt [4]{c}}+x^2} \, dx}{16 a c}+\frac {\left (\frac {3 \sqrt {c} d}{\sqrt {a}}+e\right ) \int \frac {1}{\frac {\sqrt {a}}{\sqrt {c}}+\frac {\sqrt {2} \sqrt [4]{a} x}{\sqrt [4]{c}}+x^2} \, dx}{16 a c}-\frac {\left (3 \sqrt {c} d-\sqrt {a} e\right ) \int \frac {\frac {\sqrt {2} \sqrt [4]{a}}{\sqrt [4]{c}}+2 x}{-\frac {\sqrt {a}}{\sqrt {c}}-\frac {\sqrt {2} \sqrt [4]{a} x}{\sqrt [4]{c}}-x^2} \, dx}{16 \sqrt {2} a^{7/4} c^{3/4}}-\frac {\left (3 \sqrt {c} d-\sqrt {a} e\right ) \int \frac {\frac {\sqrt {2} \sqrt [4]{a}}{\sqrt [4]{c}}-2 x}{-\frac {\sqrt {a}}{\sqrt {c}}+\frac {\sqrt {2} \sqrt [4]{a} x}{\sqrt [4]{c}}-x^2} \, dx}{16 \sqrt {2} a^{7/4} c^{3/4}} \\ & = \frac {x \left (d+e x^2\right )}{4 a \left (a+c x^4\right )}-\frac {\left (3 \sqrt {c} d-\sqrt {a} e\right ) \log \left (\sqrt {a}-\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{16 \sqrt {2} a^{7/4} c^{3/4}}+\frac {\left (3 \sqrt {c} d-\sqrt {a} e\right ) \log \left (\sqrt {a}+\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{16 \sqrt {2} a^{7/4} c^{3/4}}+\frac {\left (3 \sqrt {c} d+\sqrt {a} e\right ) \text {Subst}\left (\int \frac {1}{-1-x^2} \, dx,x,1-\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt {2} a^{7/4} c^{3/4}}-\frac {\left (3 \sqrt {c} d+\sqrt {a} e\right ) \text {Subst}\left (\int \frac {1}{-1-x^2} \, dx,x,1+\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt {2} a^{7/4} c^{3/4}} \\ & = \frac {x \left (d+e x^2\right )}{4 a \left (a+c x^4\right )}-\frac {\left (3 \sqrt {c} d+\sqrt {a} e\right ) \tan ^{-1}\left (1-\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt {2} a^{7/4} c^{3/4}}+\frac {\left (3 \sqrt {c} d+\sqrt {a} e\right ) \tan ^{-1}\left (1+\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{8 \sqrt {2} a^{7/4} c^{3/4}}-\frac {\left (3 \sqrt {c} d-\sqrt {a} e\right ) \log \left (\sqrt {a}-\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{16 \sqrt {2} a^{7/4} c^{3/4}}+\frac {\left (3 \sqrt {c} d-\sqrt {a} e\right ) \log \left (\sqrt {a}+\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{16 \sqrt {2} a^{7/4} c^{3/4}} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.20 (sec) , antiderivative size = 267, normalized size of antiderivative = 0.97 \[ \int \frac {d+e x^2}{\left (a+c x^4\right )^2} \, dx=\frac {\frac {8 a x \left (d+e x^2\right )}{a+c x^4}-\frac {2 \sqrt {2} \sqrt [4]{a} \left (3 \sqrt {c} d+\sqrt {a} e\right ) \arctan \left (1-\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{c^{3/4}}+\frac {2 \sqrt {2} \sqrt [4]{a} \left (3 \sqrt {c} d+\sqrt {a} e\right ) \arctan \left (1+\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{c^{3/4}}+\frac {\sqrt {2} \left (-3 \sqrt [4]{a} \sqrt {c} d+a^{3/4} e\right ) \log \left (\sqrt {a}-\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{c^{3/4}}+\frac {\sqrt {2} \left (3 \sqrt [4]{a} \sqrt {c} d-a^{3/4} e\right ) \log \left (\sqrt {a}+\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{c^{3/4}}}{32 a^2} \]

[In]

Integrate[(d + e*x^2)/(a + c*x^4)^2,x]

[Out]

((8*a*x*(d + e*x^2))/(a + c*x^4) - (2*Sqrt[2]*a^(1/4)*(3*Sqrt[c]*d + Sqrt[a]*e)*ArcTan[1 - (Sqrt[2]*c^(1/4)*x)
/a^(1/4)])/c^(3/4) + (2*Sqrt[2]*a^(1/4)*(3*Sqrt[c]*d + Sqrt[a]*e)*ArcTan[1 + (Sqrt[2]*c^(1/4)*x)/a^(1/4)])/c^(
3/4) + (Sqrt[2]*(-3*a^(1/4)*Sqrt[c]*d + a^(3/4)*e)*Log[Sqrt[a] - Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2])/c^(
3/4) + (Sqrt[2]*(3*a^(1/4)*Sqrt[c]*d - a^(3/4)*e)*Log[Sqrt[a] + Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2])/c^(3
/4))/(32*a^2)

Maple [C] (verified)

Result contains higher order function than in optimal. Order 9 vs. order 3.

Time = 0.24 (sec) , antiderivative size = 67, normalized size of antiderivative = 0.24

method result size
risch \(\frac {\frac {e \,x^{3}}{4 a}+\frac {d x}{4 a}}{c \,x^{4}+a}+\frac {\munderset {\textit {\_R} =\operatorname {RootOf}\left (c \,\textit {\_Z}^{4}+a \right )}{\sum }\frac {\left (\textit {\_R}^{2} e +3 d \right ) \ln \left (x -\textit {\_R} \right )}{\textit {\_R}^{3}}}{16 a c}\) \(67\)
default \(d \left (\frac {x}{4 a \left (c \,x^{4}+a \right )}+\frac {3 \left (\frac {a}{c}\right )^{\frac {1}{4}} \sqrt {2}\, \left (\ln \left (\frac {x^{2}+\left (\frac {a}{c}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {a}{c}}}{x^{2}-\left (\frac {a}{c}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {a}{c}}}\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {a}{c}\right )^{\frac {1}{4}}}+1\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {a}{c}\right )^{\frac {1}{4}}}-1\right )\right )}{32 a^{2}}\right )+e \left (\frac {x^{3}}{4 a \left (c \,x^{4}+a \right )}+\frac {\sqrt {2}\, \left (\ln \left (\frac {x^{2}-\left (\frac {a}{c}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {a}{c}}}{x^{2}+\left (\frac {a}{c}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {a}{c}}}\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {a}{c}\right )^{\frac {1}{4}}}+1\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {a}{c}\right )^{\frac {1}{4}}}-1\right )\right )}{32 a c \left (\frac {a}{c}\right )^{\frac {1}{4}}}\right )\) \(245\)

[In]

int((e*x^2+d)/(c*x^4+a)^2,x,method=_RETURNVERBOSE)

[Out]

(1/4*e/a*x^3+1/4*d/a*x)/(c*x^4+a)+1/16/a/c*sum((_R^2*e+3*d)/_R^3*ln(x-_R),_R=RootOf(_Z^4*c+a))

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 873 vs. \(2 (192) = 384\).

Time = 0.29 (sec) , antiderivative size = 873, normalized size of antiderivative = 3.17 \[ \int \frac {d+e x^2}{\left (a+c x^4\right )^2} \, dx=\frac {4 \, e x^{3} - {\left (a c x^{4} + a^{2}\right )} \sqrt {-\frac {a^{3} c \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} + 6 \, d e}{a^{3} c}} \log \left (-{\left (81 \, c^{2} d^{4} - a^{2} e^{4}\right )} x + {\left (a^{6} c^{2} e \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} + 27 \, a^{2} c^{2} d^{3} - 3 \, a^{3} c d e^{2}\right )} \sqrt {-\frac {a^{3} c \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} + 6 \, d e}{a^{3} c}}\right ) + {\left (a c x^{4} + a^{2}\right )} \sqrt {-\frac {a^{3} c \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} + 6 \, d e}{a^{3} c}} \log \left (-{\left (81 \, c^{2} d^{4} - a^{2} e^{4}\right )} x - {\left (a^{6} c^{2} e \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} + 27 \, a^{2} c^{2} d^{3} - 3 \, a^{3} c d e^{2}\right )} \sqrt {-\frac {a^{3} c \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} + 6 \, d e}{a^{3} c}}\right ) + {\left (a c x^{4} + a^{2}\right )} \sqrt {\frac {a^{3} c \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} - 6 \, d e}{a^{3} c}} \log \left (-{\left (81 \, c^{2} d^{4} - a^{2} e^{4}\right )} x + {\left (a^{6} c^{2} e \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} - 27 \, a^{2} c^{2} d^{3} + 3 \, a^{3} c d e^{2}\right )} \sqrt {\frac {a^{3} c \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} - 6 \, d e}{a^{3} c}}\right ) - {\left (a c x^{4} + a^{2}\right )} \sqrt {\frac {a^{3} c \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} - 6 \, d e}{a^{3} c}} \log \left (-{\left (81 \, c^{2} d^{4} - a^{2} e^{4}\right )} x - {\left (a^{6} c^{2} e \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} - 27 \, a^{2} c^{2} d^{3} + 3 \, a^{3} c d e^{2}\right )} \sqrt {\frac {a^{3} c \sqrt {-\frac {81 \, c^{2} d^{4} - 18 \, a c d^{2} e^{2} + a^{2} e^{4}}{a^{7} c^{3}}} - 6 \, d e}{a^{3} c}}\right ) + 4 \, d x}{16 \, {\left (a c x^{4} + a^{2}\right )}} \]

[In]

integrate((e*x^2+d)/(c*x^4+a)^2,x, algorithm="fricas")

[Out]

1/16*(4*e*x^3 - (a*c*x^4 + a^2)*sqrt(-(a^3*c*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3)) + 6*d*e)
/(a^3*c))*log(-(81*c^2*d^4 - a^2*e^4)*x + (a^6*c^2*e*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3))
+ 27*a^2*c^2*d^3 - 3*a^3*c*d*e^2)*sqrt(-(a^3*c*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3)) + 6*d*
e)/(a^3*c))) + (a*c*x^4 + a^2)*sqrt(-(a^3*c*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3)) + 6*d*e)/
(a^3*c))*log(-(81*c^2*d^4 - a^2*e^4)*x - (a^6*c^2*e*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3)) +
 27*a^2*c^2*d^3 - 3*a^3*c*d*e^2)*sqrt(-(a^3*c*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3)) + 6*d*e
)/(a^3*c))) + (a*c*x^4 + a^2)*sqrt((a^3*c*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3)) - 6*d*e)/(a
^3*c))*log(-(81*c^2*d^4 - a^2*e^4)*x + (a^6*c^2*e*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3)) - 2
7*a^2*c^2*d^3 + 3*a^3*c*d*e^2)*sqrt((a^3*c*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3)) - 6*d*e)/(
a^3*c))) - (a*c*x^4 + a^2)*sqrt((a^3*c*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3)) - 6*d*e)/(a^3*
c))*log(-(81*c^2*d^4 - a^2*e^4)*x - (a^6*c^2*e*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3)) - 27*a
^2*c^2*d^3 + 3*a^3*c*d*e^2)*sqrt((a^3*c*sqrt(-(81*c^2*d^4 - 18*a*c*d^2*e^2 + a^2*e^4)/(a^7*c^3)) - 6*d*e)/(a^3
*c))) + 4*d*x)/(a*c*x^4 + a^2)

Sympy [A] (verification not implemented)

Time = 0.45 (sec) , antiderivative size = 136, normalized size of antiderivative = 0.49 \[ \int \frac {d+e x^2}{\left (a+c x^4\right )^2} \, dx=\operatorname {RootSum} {\left (65536 t^{4} a^{7} c^{3} + 3072 t^{2} a^{4} c^{2} d e + a^{2} e^{4} + 18 a c d^{2} e^{2} + 81 c^{2} d^{4}, \left ( t \mapsto t \log {\left (x + \frac {4096 t^{3} a^{6} c^{2} e + 144 t a^{3} c d e^{2} - 432 t a^{2} c^{2} d^{3}}{a^{2} e^{4} - 81 c^{2} d^{4}} \right )} \right )\right )} + \frac {d x + e x^{3}}{4 a^{2} + 4 a c x^{4}} \]

[In]

integrate((e*x**2+d)/(c*x**4+a)**2,x)

[Out]

RootSum(65536*_t**4*a**7*c**3 + 3072*_t**2*a**4*c**2*d*e + a**2*e**4 + 18*a*c*d**2*e**2 + 81*c**2*d**4, Lambda
(_t, _t*log(x + (4096*_t**3*a**6*c**2*e + 144*_t*a**3*c*d*e**2 - 432*_t*a**2*c**2*d**3)/(a**2*e**4 - 81*c**2*d
**4)))) + (d*x + e*x**3)/(4*a**2 + 4*a*c*x**4)

Maxima [A] (verification not implemented)

none

Time = 0.28 (sec) , antiderivative size = 253, normalized size of antiderivative = 0.92 \[ \int \frac {d+e x^2}{\left (a+c x^4\right )^2} \, dx=\frac {e x^{3} + d x}{4 \, {\left (a c x^{4} + a^{2}\right )}} + \frac {\frac {2 \, \sqrt {2} {\left (3 \, \sqrt {c} d + \sqrt {a} e\right )} \arctan \left (\frac {\sqrt {2} {\left (2 \, \sqrt {c} x + \sqrt {2} a^{\frac {1}{4}} c^{\frac {1}{4}}\right )}}{2 \, \sqrt {\sqrt {a} \sqrt {c}}}\right )}{\sqrt {a} \sqrt {\sqrt {a} \sqrt {c}} \sqrt {c}} + \frac {2 \, \sqrt {2} {\left (3 \, \sqrt {c} d + \sqrt {a} e\right )} \arctan \left (\frac {\sqrt {2} {\left (2 \, \sqrt {c} x - \sqrt {2} a^{\frac {1}{4}} c^{\frac {1}{4}}\right )}}{2 \, \sqrt {\sqrt {a} \sqrt {c}}}\right )}{\sqrt {a} \sqrt {\sqrt {a} \sqrt {c}} \sqrt {c}} + \frac {\sqrt {2} {\left (3 \, \sqrt {c} d - \sqrt {a} e\right )} \log \left (\sqrt {c} x^{2} + \sqrt {2} a^{\frac {1}{4}} c^{\frac {1}{4}} x + \sqrt {a}\right )}{a^{\frac {3}{4}} c^{\frac {3}{4}}} - \frac {\sqrt {2} {\left (3 \, \sqrt {c} d - \sqrt {a} e\right )} \log \left (\sqrt {c} x^{2} - \sqrt {2} a^{\frac {1}{4}} c^{\frac {1}{4}} x + \sqrt {a}\right )}{a^{\frac {3}{4}} c^{\frac {3}{4}}}}{32 \, a} \]

[In]

integrate((e*x^2+d)/(c*x^4+a)^2,x, algorithm="maxima")

[Out]

1/4*(e*x^3 + d*x)/(a*c*x^4 + a^2) + 1/32*(2*sqrt(2)*(3*sqrt(c)*d + sqrt(a)*e)*arctan(1/2*sqrt(2)*(2*sqrt(c)*x
+ sqrt(2)*a^(1/4)*c^(1/4))/sqrt(sqrt(a)*sqrt(c)))/(sqrt(a)*sqrt(sqrt(a)*sqrt(c))*sqrt(c)) + 2*sqrt(2)*(3*sqrt(
c)*d + sqrt(a)*e)*arctan(1/2*sqrt(2)*(2*sqrt(c)*x - sqrt(2)*a^(1/4)*c^(1/4))/sqrt(sqrt(a)*sqrt(c)))/(sqrt(a)*s
qrt(sqrt(a)*sqrt(c))*sqrt(c)) + sqrt(2)*(3*sqrt(c)*d - sqrt(a)*e)*log(sqrt(c)*x^2 + sqrt(2)*a^(1/4)*c^(1/4)*x
+ sqrt(a))/(a^(3/4)*c^(3/4)) - sqrt(2)*(3*sqrt(c)*d - sqrt(a)*e)*log(sqrt(c)*x^2 - sqrt(2)*a^(1/4)*c^(1/4)*x +
 sqrt(a))/(a^(3/4)*c^(3/4)))/a

Giac [A] (verification not implemented)

none

Time = 0.28 (sec) , antiderivative size = 268, normalized size of antiderivative = 0.97 \[ \int \frac {d+e x^2}{\left (a+c x^4\right )^2} \, dx=\frac {e x^{3} + d x}{4 \, {\left (c x^{4} + a\right )} a} + \frac {\sqrt {2} {\left (3 \, \left (a c^{3}\right )^{\frac {1}{4}} c^{2} d + \left (a c^{3}\right )^{\frac {3}{4}} e\right )} \arctan \left (\frac {\sqrt {2} {\left (2 \, x + \sqrt {2} \left (\frac {a}{c}\right )^{\frac {1}{4}}\right )}}{2 \, \left (\frac {a}{c}\right )^{\frac {1}{4}}}\right )}{16 \, a^{2} c^{3}} + \frac {\sqrt {2} {\left (3 \, \left (a c^{3}\right )^{\frac {1}{4}} c^{2} d + \left (a c^{3}\right )^{\frac {3}{4}} e\right )} \arctan \left (\frac {\sqrt {2} {\left (2 \, x - \sqrt {2} \left (\frac {a}{c}\right )^{\frac {1}{4}}\right )}}{2 \, \left (\frac {a}{c}\right )^{\frac {1}{4}}}\right )}{16 \, a^{2} c^{3}} + \frac {\sqrt {2} {\left (3 \, \left (a c^{3}\right )^{\frac {1}{4}} c^{2} d - \left (a c^{3}\right )^{\frac {3}{4}} e\right )} \log \left (x^{2} + \sqrt {2} x \left (\frac {a}{c}\right )^{\frac {1}{4}} + \sqrt {\frac {a}{c}}\right )}{32 \, a^{2} c^{3}} - \frac {\sqrt {2} {\left (3 \, \left (a c^{3}\right )^{\frac {1}{4}} c^{2} d - \left (a c^{3}\right )^{\frac {3}{4}} e\right )} \log \left (x^{2} - \sqrt {2} x \left (\frac {a}{c}\right )^{\frac {1}{4}} + \sqrt {\frac {a}{c}}\right )}{32 \, a^{2} c^{3}} \]

[In]

integrate((e*x^2+d)/(c*x^4+a)^2,x, algorithm="giac")

[Out]

1/4*(e*x^3 + d*x)/((c*x^4 + a)*a) + 1/16*sqrt(2)*(3*(a*c^3)^(1/4)*c^2*d + (a*c^3)^(3/4)*e)*arctan(1/2*sqrt(2)*
(2*x + sqrt(2)*(a/c)^(1/4))/(a/c)^(1/4))/(a^2*c^3) + 1/16*sqrt(2)*(3*(a*c^3)^(1/4)*c^2*d + (a*c^3)^(3/4)*e)*ar
ctan(1/2*sqrt(2)*(2*x - sqrt(2)*(a/c)^(1/4))/(a/c)^(1/4))/(a^2*c^3) + 1/32*sqrt(2)*(3*(a*c^3)^(1/4)*c^2*d - (a
*c^3)^(3/4)*e)*log(x^2 + sqrt(2)*x*(a/c)^(1/4) + sqrt(a/c))/(a^2*c^3) - 1/32*sqrt(2)*(3*(a*c^3)^(1/4)*c^2*d -
(a*c^3)^(3/4)*e)*log(x^2 - sqrt(2)*x*(a/c)^(1/4) + sqrt(a/c))/(a^2*c^3)

Mupad [B] (verification not implemented)

Time = 14.34 (sec) , antiderivative size = 637, normalized size of antiderivative = 2.32 \[ \int \frac {d+e x^2}{\left (a+c x^4\right )^2} \, dx=\frac {\frac {e\,x^3}{4\,a}+\frac {d\,x}{4\,a}}{c\,x^4+a}-2\,\mathrm {atanh}\left (\frac {c^2\,e^2\,x\,\sqrt {\frac {e^2\,\sqrt {-a^7\,c^3}}{256\,a^6\,c^3}-\frac {9\,d^2\,\sqrt {-a^7\,c^3}}{256\,a^7\,c^2}-\frac {3\,d\,e}{128\,a^3\,c}}}{2\,\left (\frac {c\,e^3}{32\,a}-\frac {9\,c^2\,d^2\,e}{32\,a^2}-\frac {27\,c\,d^3\,\sqrt {-a^7\,c^3}}{32\,a^6}+\frac {3\,d\,e^2\,\sqrt {-a^7\,c^3}}{32\,a^5}\right )}-\frac {9\,c^3\,d^2\,x\,\sqrt {\frac {e^2\,\sqrt {-a^7\,c^3}}{256\,a^6\,c^3}-\frac {9\,d^2\,\sqrt {-a^7\,c^3}}{256\,a^7\,c^2}-\frac {3\,d\,e}{128\,a^3\,c}}}{2\,\left (\frac {c\,e^3}{32}-\frac {9\,c^2\,d^2\,e}{32\,a}-\frac {27\,c\,d^3\,\sqrt {-a^7\,c^3}}{32\,a^5}+\frac {3\,d\,e^2\,\sqrt {-a^7\,c^3}}{32\,a^4}\right )}\right )\,\sqrt {-\frac {9\,c\,d^2\,\sqrt {-a^7\,c^3}-a\,e^2\,\sqrt {-a^7\,c^3}+6\,a^4\,c^2\,d\,e}{256\,a^7\,c^3}}-2\,\mathrm {atanh}\left (\frac {c^2\,e^2\,x\,\sqrt {\frac {9\,d^2\,\sqrt {-a^7\,c^3}}{256\,a^7\,c^2}-\frac {3\,d\,e}{128\,a^3\,c}-\frac {e^2\,\sqrt {-a^7\,c^3}}{256\,a^6\,c^3}}}{2\,\left (\frac {c\,e^3}{32\,a}-\frac {9\,c^2\,d^2\,e}{32\,a^2}+\frac {27\,c\,d^3\,\sqrt {-a^7\,c^3}}{32\,a^6}-\frac {3\,d\,e^2\,\sqrt {-a^7\,c^3}}{32\,a^5}\right )}-\frac {9\,c^3\,d^2\,x\,\sqrt {\frac {9\,d^2\,\sqrt {-a^7\,c^3}}{256\,a^7\,c^2}-\frac {3\,d\,e}{128\,a^3\,c}-\frac {e^2\,\sqrt {-a^7\,c^3}}{256\,a^6\,c^3}}}{2\,\left (\frac {c\,e^3}{32}-\frac {9\,c^2\,d^2\,e}{32\,a}+\frac {27\,c\,d^3\,\sqrt {-a^7\,c^3}}{32\,a^5}-\frac {3\,d\,e^2\,\sqrt {-a^7\,c^3}}{32\,a^4}\right )}\right )\,\sqrt {-\frac {a\,e^2\,\sqrt {-a^7\,c^3}-9\,c\,d^2\,\sqrt {-a^7\,c^3}+6\,a^4\,c^2\,d\,e}{256\,a^7\,c^3}} \]

[In]

int((d + e*x^2)/(a + c*x^4)^2,x)

[Out]

((e*x^3)/(4*a) + (d*x)/(4*a))/(a + c*x^4) - 2*atanh((c^2*e^2*x*((e^2*(-a^7*c^3)^(1/2))/(256*a^6*c^3) - (9*d^2*
(-a^7*c^3)^(1/2))/(256*a^7*c^2) - (3*d*e)/(128*a^3*c))^(1/2))/(2*((c*e^3)/(32*a) - (9*c^2*d^2*e)/(32*a^2) - (2
7*c*d^3*(-a^7*c^3)^(1/2))/(32*a^6) + (3*d*e^2*(-a^7*c^3)^(1/2))/(32*a^5))) - (9*c^3*d^2*x*((e^2*(-a^7*c^3)^(1/
2))/(256*a^6*c^3) - (9*d^2*(-a^7*c^3)^(1/2))/(256*a^7*c^2) - (3*d*e)/(128*a^3*c))^(1/2))/(2*((c*e^3)/32 - (9*c
^2*d^2*e)/(32*a) - (27*c*d^3*(-a^7*c^3)^(1/2))/(32*a^5) + (3*d*e^2*(-a^7*c^3)^(1/2))/(32*a^4))))*(-(9*c*d^2*(-
a^7*c^3)^(1/2) - a*e^2*(-a^7*c^3)^(1/2) + 6*a^4*c^2*d*e)/(256*a^7*c^3))^(1/2) - 2*atanh((c^2*e^2*x*((9*d^2*(-a
^7*c^3)^(1/2))/(256*a^7*c^2) - (3*d*e)/(128*a^3*c) - (e^2*(-a^7*c^3)^(1/2))/(256*a^6*c^3))^(1/2))/(2*((c*e^3)/
(32*a) - (9*c^2*d^2*e)/(32*a^2) + (27*c*d^3*(-a^7*c^3)^(1/2))/(32*a^6) - (3*d*e^2*(-a^7*c^3)^(1/2))/(32*a^5)))
 - (9*c^3*d^2*x*((9*d^2*(-a^7*c^3)^(1/2))/(256*a^7*c^2) - (3*d*e)/(128*a^3*c) - (e^2*(-a^7*c^3)^(1/2))/(256*a^
6*c^3))^(1/2))/(2*((c*e^3)/32 - (9*c^2*d^2*e)/(32*a) + (27*c*d^3*(-a^7*c^3)^(1/2))/(32*a^5) - (3*d*e^2*(-a^7*c
^3)^(1/2))/(32*a^4))))*(-(a*e^2*(-a^7*c^3)^(1/2) - 9*c*d^2*(-a^7*c^3)^(1/2) + 6*a^4*c^2*d*e)/(256*a^7*c^3))^(1
/2)